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Abstract. We prove that the cosine law for spherical triangles and spherical tetrahedra 
defines intcgrable systems, both in the sense of multidimensional consistency and in the 
sense of dynamical systems. 

1. Introduction 

Nowadays, it is a well accepted fact that many fundamental notions and results of dif- 
ferential geometry are closely related with the theory of integrable systems. Moreover, the 
recent development of discrete differential geometry [l] led to uncovering and identifying 
the true roots of such relations in the most fundamental incidence theorems of projective 
geometry and its metric subgeometries. On the other hand, this development led also 
to a clarification of the notion of integrability, which is conceptually much more clear at 
the level of discrete systems. In particular, the notion of multidimensional consistency as 
integrability crystallized in the course of this development (see [l]). 

In the present paper, we offer a new evidence in favor of the thesis that the very basic 
structures of geometry are integrable in their nature. We show that the cosine laws for 
spherical simplices define integrable systems in several senses. 

Our main results follows: 

• We show that the map sending planar angles of a spherical triangle to its side 
lengths is integrable in the sense of multidimensional consistency (Sect. [3]). More- 
over, we identify this map with the symmetric reduction of the well known discrete 
Darboux system, which appeared for the first time in the context of discrete con- 
jugate nets. We also discuss the relation of this result with a geometric realization 
of tetrahedral angle structures, due to Feng Luo |jllj. 

• We show that the same map, interpreted as a dynamical system, is completely 
integrable in the Liouville- Arnold sense (Sect. [4]). This map can be considered 
as a non-rational discretization of the famous Euler top, which is a basic example 
of integrable system. Moreover, we identify the second iterate of this map as 
a birational discretization of the Euler top, introduced by Hirota and Kimura 
(5) and studied in detail in our previous publication |15|. This second iterate is 
shown to admit a nice geometric interpretation as a particularly simple and nice 
map between spherical triangles, which we call switch (the angles of the switched 
triangle are set to be equal to the sides of the original one). A closely related 
involutive map between spherical triangles was studied in a paper by Jonas [7], 
which can be thus considered as an early instance of appearance of integrable 
systems in elementary geometry. 
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• We show that the map sending dihedral angles of a spherical tetrahedron to its side 
lengths is a further example of an integrable dynamical system (Sect. [5]). For this 
map, we establish conserved quantities, which turn out to be a re-interpretation 
of the sine laws, and an invariant volume form. The latter is shown to be closely 
related, via the famous Schlafli formula, to the determinant of the Hesse matrix 
of the volume of the spherical tetrahedron as a function of dihedral angles. 

For convenience of the reader, we recall (in Sect. [2]) the standard derivation of the spherical 
cosine and sine laws via Gram matrices. 

2. Cosine law for spherical simplices 

The maps we are dealing with in the present paper have their origin in the cosine law 
for spherical simplices. We understand a spherical (n — l)-simplex as an intersection of 
n hemispheres in § n_1 in general position. Let Vi,...,v n G § n_1 be the vertices of a 
spherical (n — l)-simplex, and let v{, . . . ,v* G § n_1 be the vertices of the polar simplex, 
defined by the conditions 

(vlvj) = ifi^j, (v*, Vi )>0. 
The edge lengths £ij G (0, n) and the dihedral angles aiy G (0, it) are given by 

COS tij = (Vi ,Vj), COS OCij = - (v* , V*) . 

Each set of data {%} or {ay} defines a spherical simplex uniquely, and the spherical 
cosine law gives a transition from one set of data to another. The most direct way to 
derive the cosine law is through the so called Gram matrices 

G = (- cosay)£ j=1 , G' = (cos £ij)l j=1 . 



According to Lemma 1.2 in 10 , Gram matrices of spherical simplices are characterized 
by the following properties: they are positive definite symmetric matrices with diagonal 
entries equal to 1. Introducing the matrices V and W whose columns are Vi, resp. v*, we 
can write: 

W T W = G, V T V = G', V T W = D, 
where D is a diagonal matrix with positive entries. There follows: 

G' = DG~ X D, (1) 

from which the cosine laws follow immediately. They can be put as 

a'. 

cos £jj = — jg — , cosay = (2) 
where g^, g'^ are the cofactors of the Gram matrices G, G', respectively, see, for instance, 



18 . We now consider in more detail the particular cases n = 3, 4 of this construction. 



2.1. Spherical triangles. For a spherical triangle, is just the plane (inner) angle 
opposite to the edge £y, see Fig. [TJ 

The angles of a spherical triangle belong to the tetrahedron-shaped domain 

& = {a = (cti2, cti 3 , 0*23) G (0, it) 3 : a 12 + a 13 + a 2 s > 7T, -ay + a ik + a jk < n}. 

Here and below we use k) to denote any permutation of (1,2,3). The sides belong to 
the domain 

ST* = {£ = (£ 12 , £ 13 , £ 23 ) G (0, tt) 3 : £ l2 + £ l3 + £ 23 < 2tt, -i ld + £ lk + £ jk > 0}, 
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Figure 1. A spherical triangle. 



which can be obtained from via the reflection / : -h- ^"*, 1(a) = tt — a. The space 
of spherical triangles (up to orthogonal transformations of 1R 3 ) is parametrized by either 
of 3F or a spherical triangle exists if and only if a G resp. if and only if I e . 
The cosine laws for spherical triangles, as given in (pi), read: 

cos OLij + cos cos a jk , n . 

COS^j = J - : — , (3) 

sin an- sin 

cos in — cos cos £jk , , N 

cosaij = ^— — — J —. 4 

sm4sm^ 

We will consider formulas as defining the mutually inverse maps 

$ : & ->■ : ^* ->■ ^, 

a i-)- I, £ i — y a. 

We note that = / o $ o I. 

The cosine law can be supplemented by the sine law for spherical triangles: 

sin£i2 sin£i 3 sin£ 2 3 V" Y /rN 

-?T7' ( 5 ) 



smai2 sm«i 3 sina 2 3 7 



where 



and 



(i = detG = 1 — cos 2 «i2 — cos 2 — cos 2 a 23 — 2 cos cosai3 cos a 23 , 
</ = det G' — 1 — cos 2 f 12 — cos 2 £13 — cos 2 £23 + 2 cos £12 cos £13 cos £23, 



7 = sin CK12 sin ai3 sin 023, 7 = sin ti 2 sin «i 3 sin c 23 . 



The sine law is derived by considering principal 2x2 minors of the relation (fl|. 
Upon the change of variables 

COS 01 ij , Uij COS £ij , 
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the cosine law for spherical triangles takes the form 



Vij = r J ^=> ( 6 ) 



1 %2 jk 



We denote these maps by 



x _ Vij VikVjk /^\ 

-yl\l l -y 2 ik 



: t — > t*, 4> 1 : r* — > r, 

x i— > y x, 



and call them the algebraic forms of $, respectively. They are defined on the corre- 
sponding subsets r, r* of [—1, l] 3 , and are related by _1 = i o o z, i(x) = — x. The sine 
law in these new variables reads: 



1-3/12 1-3/13 1-3/Ia rf (7 



A2 



1 - x? 2 1 - xf 3 1 - x\\ 7 2 d' 

where 

d = 1 - ^12 - x i3 ~ x ls ~ 2x12X13X23, d! = \-y\ 2 - y% 3 - + 2y 12 yi 3 y 2 3, (9) 

and 

7 2 = (1 - 4i)(l - 40(1 - 4s), {if = (1 - 3/i 2 2 )(l - 2/i 2 3 )(l - 2/2 3 )- 

For further reference, we mention that the sine law can be represented in either of the 
following forms: 

d = (1 - 3/ 2 )(l - 4)(1 - d' = (1 - 4)(! " J&X 1 - ( 10 ) 

2.2. Spherical tetrahedra. For a spherical tetrahedron, see Fig. |2j the cosine law is 
expressed by the formula ^ with 

9ij x ij %ik%jk ■Eim^jm %km{%ik%jm •^im^jk %ij%km) 5 (H) 

Qii = 1 { X jk X jm X km) 2XjfcXj m Xfc m , (12) 

where we adopt the notation x^ = cos«jj and (i,j,k,m) stands for any permutation of 
(1,2,3,4). 

These formulae are supplemented (see [2j[9]) by the sine law, which can be put as 

sin £12 sin £34 sin^sin^ sin £14 sin £23 d 7' 
sinai2sina 3 4 sinai 3 sina24 sinc^sinc^ 7 c?'' 

and by the formulas 

COS £12 COS £34 — COS £13 cos £24 COS £13 cos £24 — COS £14 cos £23 



COS «i2 COS (^34 — COS CK13 COS Ct24 COS 013 COS CK24 — COS CH14 COS 023 



cos £14 cos £ 2 3 — cos £12 cos £34 d 7' 



COS a 14 COS «23 — COS a 12 COS «34 7 d 

where d = det G, d' = det G', and 7 = y/gugwgaagu, i = vVn022033044- 
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Figure 2. A spherical tetrahedron 

3. Spherical triangles and symmetric discrete Darboux system 

We now consider a somewhat fancy combinatorial interpretation of formulas ^ for the 
cosine law for spherical triangles. In this interpretation, we combinatorially assign the 
three inner angles of a spherical triangle and its three side lengths to the six 2-faces of 
an elementary cube, in such a manner that each angle and an opposite side are assigned 
to two opposite faces of the cube. More precisely, we assign the quantities x^ to the 
three faces of a 3-dimensional cube parallel to the coordinate plane jk, and the quantities 
i/jk = TiXjk to the three opposite faces, see Fig. |3j Here, Tj stands for the unit shift in 
the z-th coordinate direction. 




Figure 3. A map on an elementary 3D cube with fields assigned to 2-faces. 

Thus, what we consider now is, in the terminology of [I], a 3D system with fields 
assigned to elementary squares, given by the formulas 

(T 3 X 12 , T 2 X 13 , T!X 23 ) = 4> (^12, %13, X23) , 
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or, explicitly, 



T 3 x 12 



^2X13 



T±x 23 



X12 + X i3 X 23 



y^X 



2 

13 . 

X13 + X 12 X 23 



2 ' 
2:-! 



1 - x\ 2 y/l- X 

x 23 + x 12 x 13 



2 
23 



(15) 



'12V ± X 13 

This 3D system will be called the symmetric discrete Darboux system. It can be seen as 



the special case (symmetric reduction 
given by 

TkXij 



x 



Xji of the general discrete Darboux system, 



x 



1.1 



T XikXkj 



(16) 



System (|16|) is well known in the theory of discrete integrable systems of geometric origin. 

(7.20)] as a gauge version of a 



(17) 



It seems to have appeared for the first time in [8| eq. 
simpler but less symmetric form of the discrete Darboux system 

rp _ Xij -\~ XikXkj 

J- kXij 



1 XkjXjk 

Either form of the discrete Darboux system describes the so called 



derived earlier in [3 
rotation coefficients of discrete conjugate nets / 



Z d — > MP™, i.e., nets with planar 



elementary quadrilaterals, see [Tj p. 42]. The version of rotation coefficients satisfying 
system (16) has an advantage of being defined in more local terms, see [TJ p. 95-96]. 



Xji , and 



Besides, system (16) admits, unlike system (17), the symmetric reduction x tJ 
is therefore better suited for the description of the so called symmetric conjugate nets and 
discrete Egorov nets, see 14]. The symmetric discrete Darboux system is a close relative 



of the discrete CKP system, see 17 



Thus, we have proved that the symmetric discrete Darboux system (15) admits an 



alternative interpretation in terms of spherical triangles. This interpretation leads also to 
a new proof of integrability of the symmetric discrete Darboux system. Let us recall that 
integrability of discrete 3D systems is synonymous with their 4D consistency fl] . The 
meaning of this notion can be described as follows. 

Consider the initial value problem with the data x^ = x^, i,j = 1,2,3,4, prescribed 
at six squares adjacent to one common vertex of the 4D cube. Then the application of 



a 3D map like (15) to the four 3D cubes adjacent to this vertex allows us to determine 
all T k Xij. At the second stage, the map is applied to the other four 3-faces of the 4D 
cube, with the result being all T m {T k Xij) computed in two different ways (since each of 
the corresponding squares is shared by two 3-faces), see Fig. [1] Now, 4D consistency of 
the map means that T m (T k Xij) = T k iT m Xij) for any permutation k, m of 1, 2, 3, 4 and 
for arbitrary initial data. 



Theorem 1. Symmetric discrete Darboux system (15) is J^D consistent 



Proof. We will show that this is just a corollary of the geometry of spherical tetrahedra. 
Consider a spherical tetrahedron in § 3 with the edge lengths lij (connecting the vertices 
i and j) and dihedral angles q;^ (at the edges (km), where k, m is the complement of 
i, j in 1,2,3,4). One can take either or 1^ as independent data parametrizing the 
spherical tetrahedron. The two sets of data are related by the cosine law for spherical 
tetrahedra, see section 12.21 What seems to be less known is that this cosine law for 
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Figure 4. Initial data and results of two-stage application of a 4D consis- 
tent map with fields on 2-faces. 

spherical tetrahedra can be derived via a two-fold application of the cosine law for the 
spherical triangles, if the combinatorics of this two-fold application is in agreement with 
Fig. |4j In other words, a spherical tetrahedron combinatorially corresponds to a 4D cube. 




Figure 5. Left: Sides and angles in the spherical triangle Lk(m). Right: 
Sides and angles in the spherical triangle A(ijk) 

To show this, we start with the dihedral angles as initial data. From these, we 
can determine all the planar angles atjj* (the latter is the notation for the planar angle 
opposite to the side (ij) in the spherical triangle A(ijk) with the vertices k, see Fig. [5] 
(right)). Indeed, draw a great 2-sphere S m C § 3 with the pole at the vertex m. It will cut 
the tetrahedron along a spherical triangle which is called Lk(m), the link of the vertex m 
(see [12] )• This triangle has, as one easily sees, the angles a^-, and ctjk and the side 
lengths a[™\ and oc^\ see Fig. [5] (left). Therefore, we find: 

/ (m) (to) (m)\ 1 / \ 

' a ik > a )k ) = ®( a ij, a ik, Ujk)- 

In the above combinatorial interpretation, the dihedral angles are assigned to six 2- 
faces of the 4D cube, adjacent to one of its vertices. The four 3-faces (3D cubes) adjacent 
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to this vertex, correspond combinatorially to the four spherical triangles Lk(m). The result 
of application of our 3D system to these four 3D cubes can be written as a\^=T k aij. 

At the second stage, we apply the cosine law to the four faces A(i, j, k) of the tetra- 
hedron, see Fig. [5] (right). In the above combinatorial interpretation, these four spherical 
triangles correspond to the four 3-faces of the cube adjacent to its opposite vertex. In 
particular, the side lengths £y are assigned to the six 2-faces of the 4D cube, adjacent to 
this opposite vertex. We find: 

(%,^,^) = $(4 ) ,«;?,««)- 

This can be written as £ij = TkOc^ = Tk(T m aij). Clearly, each edge is shared by two 
2-faces, so we get two answers for each i^. Of course, they must coincide as definite 
geometric quantities in a given spherical tetrahedron. □ 

The proof of Theorem [T] given above is completely free of computations and it is based 
on the existence of a spherical tetrahedron with given dihedral angles. We mention that 



very similar arguments were used in 11 in the opposite direction: the coincidence of the 
two answers for each £{j was demonstrated there by a direct computation and used to 
provide a geometric realization for a given angle structure. 

Remark. In the combinatorial interpretation adopted in the present section, the map 
cf) : a t— > I is interpreted as solving a spherical triangle by its three planar angles, and 
combinatorially corresponds to one (directed) long diagonal of the cube. The inverse map 
(f)^ 1 : i i— > a would be interpreted as solving a spherical triangle by its three sides and 
would combinatorially correspond to the opposite orientation of this same long diagonal. 
It is easy to realize that other long diagonals of the cube would correspond to solving a 
spherical triangle by two sides and a planar angle enclosed by them, or by a side and two 
adjacent planar angles. The corresponding maps could be called companion maps to 0. It 
is well known that companion maps for a multidimensionally consistent map with fields 
on 2-faces satisfy the so called functional tetrahedron equation, see, e.g., [1, Sect. 6.16]. 
We come to a conclusion that the cosine law for spherical triangles delivers an example of 
solutions of the functional tetrahedron equation. According to a private communication 
by R. Kashaev, he was aware of this example since many years. 

4. Spherical triangles and discrete time Euler top 

In the present section, it will be more convenient to change notation for the edges 
and angles of a spherical triangle to (^1,^2,^3) for (^23,^13,^12), resp. to (ati, a 2 , cn 3 ) for 
(£23, «i3, CK12). We study a dynamical system generated by iterations of the map (ft from 
(pi), resp. its inverse map from ([7]), which now take the form 

; y. = ^ ■ X r r l- _ -l . x . _ Vi~ VjVh ,1^ 



Due to homogeneity, one can insert a small parameter e into the map by a scaling 
transformation h-> EXi, Hi ei/i leading to 
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Geometrically, small values of e in (19) correspond to small values of x,i = cosojj in (18), 
i.e., to spherical triangles with all three angles (and all three sides) being close ton/2. The 
form (19) of the map <fi makes it obvious that it can be considered as a time discretization 



of the famous Euler top, described by the following system of differential equations: 



XjXfc. 



(20) 



The latter system can be brought to the form Xi = CiXjXk, involving arbitrary real param- 
eters Q, by a simple scaling transformation x^ i-> where, however, some of 6j may be 
imaginary, depending on the signs of q. 



It is well-known that system (|20|) is integrable in the Liouville- Arnold sense and can be 

In particular, it is bi-Hamiltonian 



solved in terms of elliptic functions, see for instance 16 
and admits integrals of motion 



xj — Xj, two of which are functionally independent. 



We will show that discretization (19) of the Euler top inherits integrability. 
Theorem 2. The map <fi admits integrals of motion 

1-xl 



E, 



x 



2 ■ 



(21) 



two of which are functionally independent. It has an invariant volume form given by 

dxi A dxi A dx3 



<p(x) 



where <p{x) is any of the functions (1 — x^)(l — xl) or (1 — x\) 2 (the quotient of any two 
such functions is an integral of motion). 

Proof. The first statement (about integrals of motion) follows immediately from the sine 
law As for the second statement (about an invariant volume form), a direct compu- 
tation leads to the following expressions: 

dm = i 

dxi 



XT 



and 



8xa 



Xfc XiXj 



1 x])\\ 



XT 




The second factor in the last formula can be interpreted as a conjugation of the Jacobi 
matrix dy / dx by the diagonal matrix with entries y/L — xf, and is therefore irrelevant for 
the computation of its determinant, which is thus equal to 



det 



dy_ 

dx 



z?)(l -*!)(!-*§) 



i y-3 yi 
V3 1 yi 



d' 



(i - xD(i - xD(i - xiy 



According to the sine law (|8j), the latter expression is equal to 

'dy\ [i-yW-vD 



det 



dx 



[l-x%l-xt 



2\ ' 



which finishes the proof. 



□ 
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Corollary 1. The map is Poisson with respect to Poisson brackets (bi-vector fields) 
obtained by contraction of the tri-vector field (p(x)d\ A 82 A 83 with any of the one-forms 
dEij(x). Explicitly, a family of compatible invariant Poisson brackets is given by 

{xi,Xj} = CjXfe(l — x|) — CjXf.(l — x 2 ), 

where C\, C2, C3 are arbitrary constants. Thus, the map <fi is bi-Hamiltonian and is 
completely integrable in the Liouville- Arnold sense. 

A somewhat unpleasant feature of the map cf) is its non-rational character, due to the 
square roots in the denominators. A striking and unexpected fact is that the second 
iterate of is a rational (actually, a birational) map, which appeared previously in the 
literature on a different occasion. 

Theorem 3. On the set r fl </> _1 (t), where the second iterate of (p is defined, this second 
iterate <fi 2 : x t— > x is a birational map given by 

Xi + 2xjx k + 

1 — (x 2 + x\ + X 2 ) — 2X1X2X3 

These functions are the (unique) solutions of the following system of equations: 

which constitute the Hirota-Kimura discretization of the Euler top 



Xi ~ ~ 1^2 1 ^2 1 ^,2\ wzrzrz. • y > 



Proof. We let 0(x) = y, 4>{y) = x. According to formulas (18), we have 

Xi = y l -y J yk ~ = yi + ym ^ 



Using (|24|), we find, on one hand, 



Xj, Xi 



Zyjyk 



and, on the other hand, 

+ ykydiyk - yiVj) + - ykyi){yk + y»yj) 2 %i/fc 



1 - yf) x n - y'V^l a/i - y]V^yl 



3 

L2 



This proves equations (23) for . These equations are linear with respect to xi,X2,X3, 
and thus they can be solved to give the explicit birational map (22). □ 

Actually, the integrability attributes of the map (f) given in Theorem [2] and Corollary [T] 
were found previously for the rational map 2 : integrals in [5] and the invariant volume 
form along with the bi-Hamiltonian structure in [IB]. We refer also to [6 13 [14] for our 
recent study on Hirota-Kimura-type discretizations. 

While the map <3> relates the data within one spherical triangle, a natural interpretation 
of the map $ 2 is given by the following map between spherical triangles. 

Definition 1. Let A be a spherical triangle with the angles a.^ and the sides £j. The 
switched triangle A = 5(A) is the spherical triangle with inner angles 
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(whenever such a triangle exists). The switched triangle is defined up to a transformation 
from 0(3). 

Clearly, S(A) exists if and only if £ e ^H$(^), i.e., if and only if a G $~ 1 ( e ^)n^. In 
this case, we have £ = $ 2 (a). Thus, the quantities Xi = cos^ are related to Xi = cosaj by 
the Hirota-Kimura discrete time Euler top equations (22) or (23). By abusing notations 
slightly, we write (a, £) = S(a,£). By induction, we have: 

S n (a,£) = 

Therefore, S n (A) is defined iffae^n n . . . n $~ n (^). For any n, the latter 

domain contains an open neighborhood of the point (tt/2, 7r/2, 7r/2). 

When talking about the switched triangles to our colleagues, we have been often asked: 
"But isn't this just the polar triangle?" The answer is: "Of course not, since the switch 
generates a non-trivial dynamical system, while the polarity is an involution" . Neverthe- 
less, it turns out that these two transformations of spherical triangles are closely related. 
We denote the angles and the sides of the polar triangle A* = P(A) by (a*, £*) = P(a, £). 
They are given by 

a* = 7i - £ h £* = 7i - an. (26) 

Clearly, in the latter interpretation, P maps ST x , c 7* to itself. Comparing formulas (25) 
and (26), we are led to the following definition. 

Definition 2. Let A be a spherical triangle with the angles cii and the sides £i. The side 
flip A = F(A) is the spherical triangle with sides 

% = (27) 

(whenever such a triangle exists). The side flip F(A) is defined up to a transformation 
from 0(3). The map F is an involution (when considered on equivalence classes modulo 
0(3)). 

An immediate consequence of the above definition is the following claim. 

Theorem 4. The switch of spherical triangles, considered as a map on equivalence classes 
modulo 0(3), can be represented as a composition of two involutions, namely as the side 
flip followed by the polarity transformation: S = P o F. 

The angles a?j of the side flipped triangle A = F(A) are given by a = n — $ 2 («). The 
quantities Xi = cos oti and x« = cos satisfy 

or, explicitly, 

^ X{ ~\~ 2,XjXk ~\- Xi( Xj -|- Xj -\- Xfc) , , 

1 — (x\ + x\ + x\) — 2xiX 2 Xz 



The map / : x h- >■ x defined by ( 29 ) is a birational involution. 

It is very easy to find integrals of this map. Indeed, from the defining relation (27) 
yielding cos£j = — cosfj and the cosine law ^ there follows: 

COS OLi + COS Oj COS Oik COS CKj + COS Oij COS OJfc 

sin q? 7 - sin sin a;,- sin a;*; 
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As a consequence, one gets three functionally independent conserved quantities of the 



map (29): 



2 i (Xj + XjXk) . 2 n d 



cos L = 7T— ttt sm 



4)(l-4) 1 (l-x 2 )(l-x 2 



with the quantity d defined by Q. Note that the quotients = sin 2 £j/ sin 2 £j depend 
only on x 2 rather than on x^ and are therefore integrals of motion for the map <fi 2 = i o f 



(these are nothing but functions (21)). 

Map (29) appeared in an old paper |7j by H. Jonas. Actually, he considered a closely 
related angle flip A i— > A of spherical triangles, defined by the relation cJj = it — ai, and 
proved that the cosines of the side lengths of the spherical triangles A, A are related by 

X i ~\~ X i ~\~ XjXfc ~\~ X jXfo — . ^30^ 



The latter map is conjugated to (28) via i : x >->■ —x. Jonas' results include integrals of 



the map (30) and its solution in terms of elliptic functions. 



5. Spherical tetrahedra and discretization of two linearly coupled 

euler tops 

We now consider the algebraic form of the cosine law for spherical tetrahedra as a 
dynamical system. Recall that it is given by 

va = -^=, (31) 



with the expressions for the polynomials gij given in (11-12). This map, which we denote 
by if), can be considered as a time discretization of the following system of ordinary 
differential equations: 

Xij XikXjk ~t~ Xi m Xj m . (32) 

Here, recall, (i,j,k,m) stands for an arbitrary permutation of (1,2,3,4). This six- 
dimensional system turns out to consist of two linearly coupled copies of the Euler 
top. Namely, each group of three variables p^ = x^- + Xk m and = x^ — x^ m with 
ij = 12, 13, 23, satisfies its own copy of the Euler top equations: 

Pij = PikPjk, qij = qikqjk- (33) 



As a consequence, equations (32) are integrable in the Liouville- Arnold sense, with four 
independent integrals of motion which can be chosen as 

Pl2 - Pl3 = 0&12 + ^34) 2 - (^13 + 2>m) 2 , p\ 3 - p\ 3 = (x 13 + X 24 ) 2 ~ (x 23 + X 14 ) 2 , 
<?12 - 9l3 = (^12 - X u) 2 - («13 - ^24) 2 , ql 3 ~ <?2 3 = (^13 - X 24 ) 2 ~ (^23 ~ X lA f \ 

or, alternatively, as 

2,2 2 2 2,2 2 2 

x 12 * x 34 x 13 x 24> x 13 ' x 24 x 23 x 14> 

^12^34 — Xl3 x 24, %13 x 24: — ^23^14- 

Theorem 5. The map if> admits four independent integrals of motion 



{1 - x 2 3 ){l - x 2 4 ) (l-x 2 3 )(l 



x 2 4 ) 



(i - x 2 2 )(i - x 2 u y (i-x 2 2 )(i-x 2 34 ) 

^12^34 — 2Jl3^24 X\ 3 X 2 A ~ ^23^14 



V/(1-X? 2 )(l-X| 4 )' V(l-^2)(1-^)' 
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It has an invariant volume form given by 

dxi2 A dxi 3 A dx 2 3 A dxu A dx 2 4 A dx^ 



<p(x) 

where (p(x) is any of the functions (1 — x 2 ) 5 / 2 (l — 



x 



km) 



\5/2 



Proof. The claim about integrals of motion follows directly from the sine laws (13-14). 
Let us turn to the proof of the claim about the invariant volume form. Our aim is to 
compute the determinant of the Jacobi matrix (dyij/dxi>j>). We use the following ordering 



of the index pairs: 12, 13, 
the following expressions: 



23, 14, 24, 34. The entries of the Jacobi matrix are given by 



dyjj 

dyjj 
dx ik 

dy^ 
dx 



x 



km 



\J 9ii9jj 

1 - T 2 

x ^km 

\J 9ii9jj 



Vjk 



9kk 
9 33 ' 



X 



km 



V 9ii9jj 



km ij 
"km 



9kk9r, 



9ii9jj 



(34) 
(35) 
(36) 



where 



1 1 

Pkrr 



VikVjm VimVjk yijUikVim VijU jkV 'jm 



Indeed, the first two expressions, ( 34-35]), are easily obtained by a direct computation, 



while the third one, (36), is a result of |12] 



The last factors in the above expressions can 

by a diagonal matrix 



be interpreted as a conjugation of the Jacobi matrix (dyij/dxiif) 

j. A 1 933 1 933 1 9 a j9u 1 9339 'u 
diag (1,4/ — ,4/ — ,4/ — ,4/ — , 4/ 

V V 922 V 911 V 922 V 9ll V 011022, 

and are therefore inessential for the computation of the determinant, 
following result: 



We come to the 



det 



dy_ 
dx 



n 



k<m 



,1 X km) 



(gil922g33944) 3/2 



1 

2/23 
2/13 
2/24 
2/14 
Pll 



2/23 
1 

2/12 
2/34 
Pll 
2/14 



2/13 
2/12 

1 

pll 

2/34 
2/24 



2/24 
2/34 

Pit 
1 

2/12 
2/13 



2/14 

Pit 

2/34 
2/12 

1 

2/23 



Pit 
2/14 
2/24 
2/13 
2/23 
1 



A direct computation with a symbolic manipulation system like Maple leads to the fol- 
lowing factorized expression for the determinant on the right-hand side: 

1 



Therefore, 



2/23 
2/13 
2/24 
2/14 
Pll 



det 



2/23 

1 

2/12 

2/34 
P24 
2/14 



dy_ 
dx 



2/13 

2/12 
1 

Pll 

2/34 
2/24 



2/24 
2/34 

Pit 
1 

2/12 
2/13 



n 



k<m 



2/14 

p£ 

2/34 
2/12 
1 

2/23 

1 — rr 



P?2 
2/14 
2/24 
2/13 
2/23 
1 



2/110222/33044 

rwi-24) 



Am/ 



n 



;i - y%>> 



(011022033044) 372 * 
(011022033044) 3/2 7'' 
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At the next step, we use the following result (proved below): 



2 \h<fl-vlY 



{9ll922g33944) 3/2 Y\ k<m {^ ~ X 2 k J 2 ' 



(37) 



As a consequence, we find with the help of (J 1 3|) 

dx 



det 



hj , 



d! 



n 



k<m\ 



X km) 7' 



d' 

i 



and one more reference to relations (13) proves the claim about the invariant volume 
form. 

To finish the proof we need to prove formula (37). Comparing expression q with @, 
we see that gu is the analogue of the quantity d for the spherical triangle with the angles 
otjk, Oij mi akm, i-e., for the triangle Lk(i). Then, from (10), we obtain: 

1/2 (i) . 

= sm a)l sm a jm sm a km . 
Multiplying all three such expressions for Lk(i), we find: 

3/2 . (i) . (i) ■ (i) ■ 2 -2 -2 

= smoA^ smavsma^sm ajfcsm a jm sxn a km . 

Similarly, g' u is the analogue of the quantity d' for the spherical triangle with the sides 
£jk, ^jmi ^km, i-e., for the triangle A(j, k,m). From (10), we obtain: 

(9'u) 1/2 = s m «5r } sin^ m sin4 m . 



Multiplying all three such expressions for A(j, k,m), we find: 
As a consequence, 



„( m ) 



(*) 



sm a , sm a:-™ sm a 



sm 



jk sm 



Sill tfc m . 



(fi'llfi'225 , 335 , 44 



,3/2 



sm a 



(0 
jk 



nn 

v«=l j<fe 

nn-»s? 

\rn=l j<k 



j<k 




sm 4 £ jk . 



Indeed, each dihedral angle oijk participates as a planar angle of two triangles Lk(i) and 
Lk(m), and each side £j k is shared by two triangles A(j, k,m) and A(i, j, k). This proves 
(|37|). □ 



Remark. According to the famous Schlafli formula (see, e.g., 12,18]), if V = V(a) 
is the volume of the spherical tetrahedron with the dihedral angles Oy, then dV = 
^kmdctij. Thus, our proof of Theorem [5] actually leads to an evaluation of the 
determinant of the Hesse matrix of the volume function, det(d 2 V / daijdai>ji) . We are not 
aware of such an evaluation in the literature. 



Remark. The reduction of the system of ordinary differential equations (32) to two 
copies of the Euler top (33) yields that the original system admits two invariant 3-forms, 
dpi2Adpi^Adp23 and A/12 A<igi 3 A<ig23- We conjecture that the discretization ip also admits 
an invariant 3-form. This property is stronger than existence of an invariant volume form 
established in Theorem [5] If true, this 3-form (or, better, the corresponding 3- vector field) 
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would yield, by contraction with gradients of integrals of motion, a family of compatible 
Poisson structures for ip. 
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